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ABSTRACT 

> 

, Whitehead's 1922 theory of gravitation continues to attract the attention of philoso- 

\ phers, despite evidence presented in 1971 that it violates experiment. We demon- 

strate that the theory strongly fails five quite different experimental tests, and conclude 
that, notwithstanding its meritorious philosophical underpinnings. Whitehead's theory 
is truly dead. 



1. Introduction and summary 



^ ■ In 1922, the distinguished mathematician and philosopher Alfred North Whitehead (1861- 

^ . 1947) , then in his 60th year, published a relativistic theory of gravity with the property, which it 

shares with Einstein's theory, of containing no arbitrary parameters. Furthermore, when suitably 
interpreted, it yields the same predictions as General Relativity (GR), not only for the three classic 
tests of light bending, gravitational redshift and the precession of the perihelion of Mercury, but 
also for the Shapiro time delay effect (Shapiro 1964), recently confirmed to one part in 10^ (Bertotti 
et al. 2003). 

The reason for this coincidence was realized early on by Eddington (1922). In the case of 
vanishing cosmological constant the Schwarzschild solution is not only an exact solution of Einstein's 
theory, it is an exact solution of Whitehead's theory as well. Thus it gives the same predictions 
for the parametrized post-Newtonian (PPN) parameters 7 = (3 = 1. Eddington's remark nicely 
explained an observation of Temple (1924) that the predictions of the precession of the perihelion 
for the two theories agree exactly, and gave rise to the (incorrect) idea that it is indistinguishable 
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from GR, a point refuted by Harvey (1964) by the observation that Birkhoff's theorem fails for 
Whitehead's theory: the field outside a spherically symmetric source is not just given by the 
Schwarzschild solution but in general contains an additional constant of integration which is in 
principle measurable. 

In fact an even stronger statement can be made. This remarkable correspondence of exact 
solutions extends to the Kerr solution (Russell & Wasserman 1987) and thus to the corresponding 
Lense-Thirring or frame dragging effects (Rayner 1955a). Thus experiments such as that involving 
the LAGEOS satellites (Ciufolini & Pavhs 2004) which have verified the effect at the 10-15% level 
and the ongoing NASA-Stanford Gravity Probe B superconducting gyroscope experiment, which 
aspires to an accuracy of 1%, cannot distinguish Whitehead's from Einstein's theory on the basis 
of frame dragging (we will see below that LAGEOS actually tests Whitehead because of the failure 
Birkhoff's theorem). 

The mathematical explanation for this striking, but accidental, coincidence is that both the 
Schwarzschild solution and the Kerr solution may be cast in Kerr-Schild form (Kerr &; Schild 1965). 
That is, coordinates exist for which 

dlJ-v ~ "I" ^jx^vi (1) 

where 

r/'^V- = 5''V<^ = 0- (2) 
and is tangent to a null geodesic congruence, 

^[Mii/^^^A] = ) (3) 

where l'^ is obtained from Ifj^ by index raising using either the metric r]ni, or the metric gfj,i,. It 
follows (Kerr & Schild 1965) that 

hjxv — l^i^v (4) 

satisfies the linearized Einstein equations. If this can be chosen to agree with Whitehead's retarded 
solution, then his metric and that of Einstein will agree exactly. 

Thus for a single particle at rest at the origin, in spherical polar Minkowski coordinates t, r, 9, (f), 
Whitehead's metric is 

ds^ = -dt^ + dr^ + r'^{de^ + sin^ ed<j)'^) + ^{dt- drf . (5) 
On the other hand, the Schwarzschild metric, in standard Schwarzschild coordinates {T,r,9,(f)) is 



ds^ = -(1 - ^)dT^ + + r\de^ + sin2 edct>'') . (6) 



If we set 

* = T-Mln(^-1), (7) 
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the two metrics are seen to coincide, which is Eddington's observation (Eddington 1922). If we 
define u = t — r, then the coordinates u,r,9,(f) are nowadays referred to as outgoing Eddington- 
Finkelstein coordinates. Thus Whitehead's spacetime manifold is gcodcsically incomplete with 
respect to his curved metric because outgoing Eddington-Finkelstein coordinates cover only the 
lower half of the full Kruskal manifold. The surface r = 2M is the past event horizon, and 
Whitehead's particle is naked and corresponds to what is now called a White Hole (Harvey 1964), 
the time reverse of a Black Hole. The Kerr solution is also of Kerr-Schild form and is also an 
exact solution of Whitehead's metric (Russell & Wasserman 1987) when expressed in terms of 
advanced null coordinates. It corresponds therefore to a naked rotating White Hole. For strong 
fields therefore, even if a single object is considered, the two theories would be expected to make very 
different predictions. Note also that this exact correspondence between solutions of Whitehead's 
theory and solutions of Einstein's theory holds only for a special class of solutions. Not every 
solution h^u of linear theory may be cast in the form (4) such that (2) and (3) hold, i.e to be of 
Kerr-Schild form. Moreover not every solution in the Kerr-Schild class need be of the retarded form 
specified by Whitehead. Thus there is no general agreement between the predictions of Einstein 
and those of Whitehead. 

In any case, for many years the two theories were considered to be experimentally indistin- 
guishable, and this gave rise to much philosophical discussion as to whether additional criteria, for 
example aesthetic considerations or philosophical preconceptions, were needed in order to reject 
or accept one of them. This is brought out in Broad's review of Whitehead's book The Principle 
of RelatAvity (Broad 1923), and a particularly clear discussion indicative of the mood in the late 
1950's is that of Bonnor (1958). 

From today's perspective, one can say that the principal difference between Einstein and 
Whitehead is the latter's insistence on fixed a priori spatio-temporal relations, which in practice 
meant the adoption of a fixed, and in particular unobservable, background Minkowski spacetime. 
This is stated with admirable clarity by the philosopher John Bain (1998), who provides a valuable 
account of how Whitehead's ideas about relativity were embedded in his overall philosophy of 
nature (see also Tanaka (1987)). 

In fact, by the late 1960's the promise of new technology had led to a more optimistic, empirical 
viewpoint, and rival theories of gravity were carefully scrutinised both for internal consistency and 
for testable predictions additional to the three classic tests. An important milestone was Shapiro's 
time delay prediction (Shapiro 1964). An outcome of this line of research was the discrediting by 
Will (1971b) of Whitehead's theory. 

However, Whitehead's philosophical ideas continue to attract widespread attention, often un- 
der the rubric of Process Philosophy, and perhaps because of his formidable achievements in the 
foundations of mathematics and logic. He was after all co-author with Bertrand Arthur William 
Russell (1872-1970) of the epoch making Principia Mathematica. As a result, many of his followers 
have been reluctant to abandon his theory of gravity despite the growing observational evidence 
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against it. 



Will's original disproof of Whitehead's theory was based on the fact that Whitehead's theory 
predicts an anisotropy in the "locally measured" Newton's constant due to distant matter. Thus a 
mass M at a distance r from the Earth produces an effective Newton's contant 



where 9 is the angle between the Earth's radial direction and the distant gravitating body. This 
would produce anomalous Earth tides that would show up in gravimetcr experiments, yet there 
was no experiment evidence for such effects (Warburton & Goodkind 1976). As a critique of 
Will's argument, it was pointed out that the resultant Earth tides depend on the distribution of 
extra-solar system matter (Mentock 1996) whose distribution is uncertain, and so a cancellation 
might take place. However, as we shall show below, allowing for these uncertainties will not 
change the predicted effect sufficiently to invalidate Will's argument. Another attempt to avoid 
Will's argument was to change the interpretation (Rcinhardt &: Rosenblum 1974; Hyman 1989). In 
Chiang & Hamity (1975) it was shown that the re-interpretation of Reinhardt &: Rosenblum (1974) 
would not achieve this goal, and they obtained the same result for the anisotropy of Newton's 
constant (8) as did Will. These general conclusions, while accepted by Bain (1998), were rejected 
by Fowler (1974) and the latter's remarks were reiterated by Tanaka (1987). Similar reservations 
have been expressed by Russell & Wasserman (1987). 

In fact, normally in science a single incorrect prediction is regarded as sufficient grounds for 
rejecting a theory; hence the well known dictum of "Darwin's Bulldog" 

The great tragedy of Science - the slaying of a beautiful hypothesis by an ugly fact. 
Thomas H. Huxley (1825 - 1895) 

By contrast, as we are reminded by Popper (1959), the confirmation of a theory is never 
complete. The best one can do is to subject it to increasingly precise and exacting tests covering 
a wider and wider range of phenomenona and circumstances. 

It turns out that Whitehead's theory is definitely excluded by several modern experiments, 
and our aim in this article and the reason for our title is to point out that any one of them is 
sufficient for rejection. In other words judged by modern scientific and technological standards, 
Whitehead's theory, beautiful as it may seem in the eyes of many of its beholders, is truly dead. 
By contrast, Einstein's theory passes all of these tests with flying colors. 

Specifically, Whitehead's theory fails five tests, most of them by many orders of magnitude 
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1. Anisotropy in G. We have reexamined Will's 1971 derivation, incorporating a model for the 
mass distribution of the galaxy that includes a dark matter halo. The predicted effect is still 
at least 100 times larger than the experimental bound. 
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2. Nordtvedt effect and lunar laser ranging. Whitehead's theory predicts that massive, self- 
gravitating bodies violate the weak equivalence principle in that their acceleration in an 
external gravitational field depends on their gravitational binding energy (Nordtvedt effect). 
The predicted size is 400 times larger than that permitted by lunar laser ranging. 

3. Gravitational radiation reaction and the binary pulsar. The theory predicts anti-damping 
of binary orbits due to gravitational radiation reaction at a level {tj/c)"^ beyond Newtonian 
gravity, in contrast to the (v/c)^ damping effect in GR. Thus it strongly violates binary pulsar 
data by about four orders of magnitude, and with the wrong sign. 

4. Violation of Birkhoff's theorem, and LAGEOS satellites. The static, spherically symmetric 
solution of the theory for finite sized bodies has an additional contribution dependent on the 
body's size (Harvey 1964; Rayner 1954; Synge 1952). This produces an additional advance of 
the perigee of the LAGEOS II satellite, in disagreement with observations by a factor of 10. 

5. Momentum conservation and the binary pulsar. Whitehead's theory predicts an acceleration 
of the center of mass of a binary system, a violation of momentum conservation (Clark 1954). 
Precise timing of the pulsar B1913+16 in the Hulse- Taylor binary pulsar rules out this effect 
by a factor of a million. 

Any of these tests alone would have been enough to kill Whitehead's theory, so collectively they 
amount to overkill. On the other hand they illustrate both the precision and depth that modern 
technology has brought to the problem of testing gravity, and serve as a warning to any would-be 
inventor of an alternative gravity theory, or to anyone who might hope that a suitably modified 
or reinterpreted Whiteheadian theory would pass muster (Schild 1956; Hyman 1989; Reinhardt &; 
Roscnblum 1974). It is not sufficient to check the "classic tests" of light bending, perihelion advance 
of Mercury, and gravitational redshift. There is now an exhaustive battery of empirical checks that 
must be done. 

The remainder of this paper provides some technical details to support these conclusions. 
Throughout, we adopt the "canonical" version of Whitehead's theory, specified as follows. One first 
assumes the presence of a flat background metric r]^^, whose Riemann tensor vanishes everywhere. 
This background metric defines null cones for any chosen spacetime event x^, given by points x'^ 
satisfying 



(9) 



The physical metric g^u is then given by (henceforth we use units in which G = c=l) 



a 

{y-y = x^-ix-y 



c/^j,(a;") = ?7^;/-2^ma 




Wa = VnuiVaYidXa/do-) 
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da^ = rii^^dx^dx" , (10) 

where the sum is over over all particles, with rest mass rua- Indices on {ua}^ are raised and lowered 
using r)niy. The quantities {x~)^^, (dx'^/da)" are to be evaluated along the past flat null cone of the 
field point x^^. 

Following Synge (1952), we assume that matter fields couple only to the physical metric g^^. 

This makes Whitehead's theory a "metric theory" of gravity (see Will (1993) for discussion). As 
such, it automatically satisfies the Einstein Equivalence Principle (EEP), which has been verified to 
extremely high precision using laboratory Eotvos-type experiments (parts in 10^^) and gravitational 
redshift experiments (parts in 10^), among others. The background metric then has no further 
direct physical consequences, apart from its role in defining the physical metric. This point mass 
expression can be generalized to continuous fluids in a straightforward way (Will 1993). 

2. Post-Newtonian limit and gravitational radiation reaction in Whitehead's theory 

2.1. Solution of Whitehead's theory to 1.5 post-Newtonian order 

We wish to evaluate the Whitehead metric within the near-zone of a slow-motion gravitating 
system, in order to derive the equations of motion. This corresponds to field points such that 
|x| ^ A ~ R/v, where A is roughly a gravitational wavelength, R and v <^ 1 are the characteristic 
size and internal velocity of the system. Accordingly, we want to evaluate g^i, at (t, x) in terms 
of source varables evaluated at the same time t. We make the standard assumption of post- 
Newtonian theory that ^ m/r ^ e, where e is a small parameter used for bookeeping purposes. 
Our goal is to determine the metric through 1.5 post-Newtonian order, or to order e^^^ beyond 
Newtonian gravity; this involves evaluating yoo through 0(e^/^), g^j through O(e^), and gij through 
0(e^/^). This will include the usual post-Newtonian terms relevant for solar-system tests, as well 
as, it will turn out, the leading effects of gravitational radiation reaction in this theory. 

We expand the retarded position of the a-th particle by 

X- = Xa(t-|x-X-|) 

~ - Va|x -X^l -h -aa|x -X^P -h . . . , (11) 

where x,j, Vq and a^ are the position, velocity and acceleration of the a-th particle at the field-point 
time t. We can then expand the spatial component (?/„)* = (x — x~y in terms of the instantaneous 
difference z^^ = {x — XaY according to 

(yaT = 4 + ^'^'viy - leaiy' + le^/'dly^ + 0{e') , (12) 

where y = ly^ |. We also expand the retarded velocity component {v~y = [dx^^/dt)^ according to 

(VaT = vi- e'/'aiy + ^edly' + 0{e'/') . (13) 
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Note that, because the quantity (y„ is null with respect to the fiat metric, (y„ )° = y, and thus 

w- = {dt/da)-{-y + e^/^a • v^) . (14) 

The foregoing expressions can then be iterated to the required order in e to convert all expressions 
into functions of v^, aa, SLa, Za, and Za = |za|. The result is 



(15) 



where 



(ha. 



00 



+ e 



2-^ 



r ■ Za , (,Vq, • Zq 
5 1 5- 



3/2 



-aa • Za - 2Va • aa + 6 



-4 



(Va • Za) 



+ 0(e2). 



(^a)oj = - 



24 



+ 2- 



+e 



ai + 4 



3-a ■ Za^a 



vizi 



(Va • Za)^2;a , Va • Zaf^ 



+ 2 



+ 0(6=^/2), 



{ha)ij — 



+ e 



1/2 ^■^a fa 



+ 0(e). 



(16) 



Indices on spatial vectors are raised and lowered using the Cartesian metric; parentheses around 
indices denote symmetrization, while square brackets denote antisymmetrization. 



The first term in (/ia)oo can be recognized as yielding the normal Newtonian potential U, given 



by 



m„ 



X ■ 



(17) 



Note the presence of 0.5PN terms in the metric; these are terms of order e^^^ in {ha)oo, and e° in 
{ha)oj- Because of general covariance, we are free to change coordinates to manipulate the form of 
the physical metric. In particular, we can remove these 0.5PN terms, can manipulate the PN terms 
to put them into a form to make comparisons with the standard parametrized post-Newtonian 
(PPN) framework (Will 1993), and can simplify the 1.5PN terms. Even though the background 
metric rj^^ will change its form under such coordinate transformations, this will have no physical 
consequences, since only ^^i^ couples to matter. 

The following coordinate transformation kills the 0.5PN terms in the physical metric, puts the 
PN terms into the standard PPN gauge, and also kills the 1.5PN terms in {ha)ij- 



t = t-_2eLO + |e3/2;^,o + 0(e2), 
x' + eX -j - 2€^/^P + 0(e2) , 



X 



(18) 
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where commas denote partial derivatives, and where 

L° = ^ ma In Za , 
a 

P = ^ ruavi In Za , 

a 

X = Y^maZa. (19) 

a 

Note that the first term in the time transformation is the post-Newtonian analogue of the Eddington 
(1922) transformation. In carrying out the normal coordinate transformation, 

to 1.5PN order, it is also necessary to reexpress the potentials in terms of the new coordinates. For 
example, the Newtonian potential changes according to 

U^U- e(C/2 + $2 + ^w) - 2e3/2[L°f7 + LWj - T.{V),j + S(t;^L°)j] , (21) 

where all quantities on the right side are barred, and where 

S(/) = "'"•^^^'''"^ • (22) 

a ^° 

The potentials $2 and are defined below. 

A further coordinate transformation, given by 

t = t'-e'^ (4?7'L'0 + 2X[^ L^o - 2MJ ) , (23) 

where = Yla''^o-Xj{'Xa) In^a, simplifies the 1.5PN terms in goj and 500- 

The post-Newtonian part of the metric will be discussed in Sec. 2.2, while the 1.5PN part will 
be discussed in Sec. 2.5. 



2.2. PPN Parameters 
Following the coordinate transformation of Eqs. (18), the metric to PN order takes the form 

500 = -1 + 2?7 - 2C/2 - 3$i - 2^>2 + 6^ - 2^w , 

goj = -AV^-'^-W^, 

gij = 6ij{l + 2U), (24) 
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where we drop the explicit use of e, and where the potentials are given by 

a ^« a ^« g,b^g 

, rUgiVg • Za)2 ^ mgnibZg ( Zgh Zb 

>l = 2^ 3 ' = 3 I 

^» a,b^a ^« ^ 

^ ^mot^^ |y^=^ "^°'^°3^°^° , (25) 

where = z^ - z^. 

We now compare this metric with the PPN metric for point masses given in Will (1993, 
2006). Transforming from the perfect fluid version of U to the point mass version using C/fjuid = 
^point — — 37$2) and working in the universal rest frame, it is a simple matter to read off the 
PPN parameters, 

7 = 1, P = l, ^ = 1, 
ai = , a2 = , as = , 

Ci = -4, C2 = -l- (26) 

Because we are dealing with point masses rather than perfect fluids, the PPN parameters ^3 and 
^4, associated with energy density and pressure are not determined. 

The parameters 7 and /? are the same as in GR. The potential is the infamous "Whitehead" 
potential, which did not appear in earlier versions of the PPN framework (Nordtvedt 1968b; Will 
1971a; Will & Nordtvedt 1972; Will 1974). It was later seen to be a generic consequence of any 
"quasi-linear" theory of gravity (Will 1973). The original PPN framework was then extended to 
incorporate naturally this potential with its associated "Whitehead parameter", ^. In Whitehead's 
theory, ^ = 1, while in GR, ^ = 0. The parameters ccj all vanish, as they do in GR, indicating that 
the theory satisfies a kind of Lorentz invariance for gravity, and has no "preferred-frame" effects. 
This is not surprising, given that it is constructed using a flat background Minkowski metric. Here 
we ignore any coupling between local gravity and a background cosmological solution for the metric, 
which can in fact lead to non-zero a's, even with a flat background metric (see Lee et al (1976) for 
an example). The "conservation-law" parameters Q are non-zero, indicating that the theory lacks 
global conservation laws for momentum and angular momentum; in GR, the C's all vanish. In Sec. 
3 we will see that many of these values are in violation of experimental bounds. 



2.3. Static spherically symmetric metric 



For a single, static point mass M, the metric is particularly simple. Placing the mass at the 
origin of coordinates, we see that (yg)^ = y = r, {ygY = , and w~ = y = r. The metric then 
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is goo = — 1 + 2M/r, goj = —2Mx^/r'^, gij = 6ij + 2Mx^x^ /r^. The coordinate transformation (7) 
converts the metric to the Schwarzschild metric of GR (Eddington 1922). This was the basis of the 
claim made in the early years of Whitehead's theory that it satisfied all the "classic" tests. 

However, real bodies such as the Sun and Earth arc not point masses, but are finite sized 
objects made up of many masses. Working in the PN limit and assuming a spherically symmetric 
collection of masses centered at the origin, it is easy to show that, for a field point outside the body. 



U = M/r , X = Mr + I/Sr , = -^2 - Mi/Sr" 



(27) 



where the latter follows from manipulating the identity V + 2U — 3^2) = ~2X^ijU^ij, and 
where M = ma and I = Yla "i^ai^a total mass and spherical moment of inertia of the 

body. All other post-Newtonian potentials vanish. Thus, in the PPN framework, the metric for a 
finite spherically symmetric static body becomes (Harvey 1964; Synge 1952) 



500 
90j 



-1 + 2M/r - 2p{M/rf + 2^MI/3r'^ 



0, 



gij = 6^J{l + 2-fM/r). 



(28) 



Recall that 7 = /3 = ^ = lin Whitehead's theory. The pericenter advance per orbit of a test 
particle moving on a geodesic of this metric is given by 



Alo = 



P 
Gvrm 

P 



(2 + 27 - /?) + 



2^ / 



1 



3 mp'^ ' 4 



•(1 + 



3 mp- 



(29) 



where p = a(l - 



with a and e being the semi-major axis and eccentricity of the orbit, and where 



the second line is the Whitehead prediction. The size-dependent term in Aui has a negligible effect 
on the perihelion advance of Mercury, and so Whitehead's theory agrees with the data; however 
that term will have measurable consequences for ranging of the Earth-orbiting LAGEOS satellites 
(Sec. 3.4). 



2.4. Anisotropy in the locally measured G 

Although we have set the fundamental coupling constant G equal to unity by making a specific 
choice of units, it turns out that, in many alternative theories, the "localh^ mcasTircd" G may 
vary. By locally measured G we mean the output of a Cavendish- type experiment, whereby one 
measures the force between a test body and a source body separated by a chosen distance. The 
result may depend on the velocity of the laboratory relative to a preferred frame, if any of the a 
PPN parameters is non-zero, and may also depend on the presence of matter outside the laboratory. 
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In the case of Whitehead's theory, there are no preferred-frame effects, but there are "preferred 
location" effects. Specificahy (see Eq. (6.75) of Wih (1993)), the locahy measured G is given by 

Giocai = 1 + ^^ext + (l - J^) e'e^C/<p , (30) 

where /, M and R are the spherical moment of inertia, mass and radius respectively of the source 
body in the Cavendish experiment, 

with the sum extending over all masses external to the laboratory, and where e* and are unit 
vectors pointing from the source body to the test body and to the a-th external body, respectively. 
Angular brackets around the indices denote a symmetric, trace-free (STF) tensor. Equation (8) 
is the special case of Eq. (30) for a single external body, and for a point source mass (7 = 0). 
The most important effect is the anisotropy in Giocab which can lead to anomalous Earth tides in 
geophysics (for the Earth, / ~ 0.5Mi?^). Notice that only the Z = 2, or quadrupole anisotropy in 
the external matter distribution contributes. 



2.5. Gravitational radiation reaction 



We focus now on the 1.5PN terms in the metric. Combining the relevant terms from Eqs. (16) 
with the 1.5PN terms generated by the coordinate transformations (18) and (23), we obtain, 



h 



(5/2) _ 
00 ~ 



.(2) _ 



vrin 



^ka ■ z„ - 2v„ • a„ + 6 ° ° - 12- 

3 



+4L°f7 - 4L^Uj + 4E(L^) J - 4E(?;-'L°) j , 



ag • ZaZi vizi ^ ^ (Va • Zafzi _ ^ Va • ZaVi 



+ 2- 



h. 



(3/2) 



-2V +4L°Uj, 
0, 



(32) 



where the superscript (n) denotes the order of e. With these expressions and the geodesic equation, 
it is straightforward to derive the 1.5PN contributions to the equation of motion of a body in the 
presence of other bodies, 

^ - l/,(5/2) _ ;^(2) _^(2) k .33^ 

We restrict attention to a binary system, and evaluate the terms in Eq. (33) at body 1 (as usual, 
dropping contributions to potentials due to body 1 itself). We use the fact that, at body 1, 
L° = m2lnr, L° = m2X^/r'^, LP = m2V2 Inr, U = m2a2lnr — m2'y2^2 -x/r^, and so on, where now 
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= x\ — X2 and r = |xi — X2I; we also recall that X^amaOa = from conservation of momentum 
at Newtonian order. The surprising result is that, despite many cancellations, there is a residual 
acceleration at 1.5PN order, given by 

TX"^ 

a{ = Smim2—r ■ 

The acceleration for body 2 is found by interchanging m\ and m2 and letting x^ 
relative acceleration = a\ — 03 is then given by 

TX^ 

a' = 16/Ltm — ^ , (35) 

where m = mi + 1112 and ^ = mim2/m are the total and reduced mass of the system, respectively. 
This radiation reaction term does not affect the orbital angular momentum, but it does cause an 
increase in the orbital energy at the rate dE/dt = 16 ix^ mr^ /r^. 

We will see in Sec. 3.3 that this has disastrous consequences for Whitehead's theory. 



(34) 
-x^. The 



2.6. Failure of momentum conservation 

In gravitational theories that lack suitable conservation laws for total momentum of gravitating 
systems, a binary system could suffer an anomalous acceleration of its center of mass, given in the 
PPN framework by 

1 m fi 5m, e 

^CM = TTC2 + a3— — — , (36) 

2 m m (1 — e^)-V^ 

where 5m = mi — m2 and the acceleration is directed toward the pericenter of the lighter body. 
In GR (and in any theory based on an invariant action) the effect vanishes, but in Whitehead's 
theory, it does not. 

Levi-Civita (1937) once claimed that this center-of-mass effect occurred in GR, but Eddington 
k, Clark (1938) spotted his error and confirmed that it did not. Clark (1954) later showed that the 
effect did occur in Whitehead's theory, in agreement with Eq. (36). At the time, of course, there 
was no hope of detecting the effect using known binary systems. However, the binary pulsar (Sec. 
3.5) provides a particularly stringent bound on this effect. 



3. Experimental tests of Whitehead's theory 

3.1. Gravimeter tests of the anisotropy in Giocai 

If Ciocai is anisotropic because of the presence of an external mass, then there will be anomalous 
tides of the solid Earth, superimposed on the normal luni-solar tides (see Nordtvedt & Will (1972); 
Will (1993) for detailed discussion). The latter are of typical amplitude /^g/g ~ 10~^ (here g is 
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the local acceleration as measured by a gravimeter) . If the external body is the sun itself, then 
Ucxt ~ 10~^, and the G anisotropy will produce a tidal signal of comparable amplitude and of 
the same frequencies as the solar tide. It is very unlikely that Whitehead's theory would survive 
a comparison between the measured solar Earth tide and standard tidal theory with such a large 
additional amplitude. However, the bound one could achieve has never been investigated in detail, 
because a cleaner test is provided by looking at the so-called sidereal tides. 

If the external mass is that of the galaxy, then ?7ext ~ 5 x 10~^, and the direction is fixed in 
space. This produces tides at frequencies associated with the sidereal day rather than the solar 
day of the solar tide, and these can be compared with known sidebands of the coupled lunar and 

solar tides. Measurements by Warburton &; Goodkind (1976) using superconducting gravimeters 
showed no evidence of anomalies, and placed the bound on the Whitehead parameter |^| < 10"^, 
as compared with the Whitehead value of unity. This improved upon the earlier bounds of Will 
(1971b), which were based on the existing tidal literature. 

This was considered a fatal blow to the theory, but it did assume an amplitude 5 x 10~^ for the 
anisotropic part of the galactic potential. That value came from relating the solar system's orbital 
velocity in the galaxy to the potential via v'^ ~ Uext- This was criticized (Mentock 1996) because 
it concentrated the mass of the galaxy at the center, whereas we now know that the bulk of the 
mass of the galaxy is in a roughly spherical halo of stars and dark matter, substantially larger in 
size than the visible Milky Way. 

However it can be shown using a simple density model for the galaxy that the original estimate 
for the anomalous tidal amplitude holds up within a factor of two. First, we note that the "trace- 
free" tensor potential = — where X is the "superpotential" defined in Eq. (19). For 
a spherically symmetric distribution of matter, X is given by 

1 f 1 f°° 

X = rm{r) + — / Airp'r^^dr' + - / ATrp'r'{r'^ + Sr'^)dr' , (37) 

3r Jo 3 Jr 

where p is the mass density and m(r) is the mass inside radius r. Then, for spherical symmetry. 



X,<ij> = h<'^>{d^X/dr^ -r-^dX/dr) 
m{r) I{r) 



(38) 



where I(r) is the spherical moment of inertia inside radius r. For flat or monotonically decreasing 
density distributions, the second term is always smaller than the first. 

To compare with the earlier estimate we consider a specific density distribution given by 
Airp = a/r'^, for r < Vc, and Airp = a/r"^, for r > rc, where Tc is a core radius meant to represent the 
mass of the inner part of the galaxy, and a is a parameter. The 1/r^ density distribution is meant 
to model the dark matter halo, and to yield a fiat rotation curve for the outer reaches of the Milky 
Way, in rough agreement with observations. By noting that a circular orbit in a spherical potential 
satisfies, v'^/r = ar = m{r)lr^^ and considering the case r > Tc, we can fit a = ^^/(l — 2g/3), and 
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find that 



T<ij> _ 2..2^<»,> l-9 + gV5 



where g = rc/r. For the case r < & similar calculation gives 

C^exP = Iv'n^'^^ , (40) 

independent of r. Thus for v 220 km/s, we find an amplitude 2 — 3 x 10^^, fully consistent with 
the earlier estimate. Note from Eq. (38) that only the matter inside our radius has an effect on 
the anisotropy. Even though the galaxy and its halo are not strictly spherically symmetric, this is 
unlikely to alter the estimate significantly. The only way to suppress this effect is by some specific, 
fine-tuned distribution of external matter. 

The conclusion stands: Whitehead's theory violates geophysical tide measurements by about 
a factor of 500. 



3.2. Lunar laser ranging and the Nordtvedt effect 

In many alternative theories of gravity, there is a violation of the weak equivalence principle 
for massive, self-gravitating bodies. Specifically, the passive gravitational mass nip may differ from 
the inertial mass according to 

nip = mi{l - ri\Eg\/mi) , (41) 

10 2 2 1 

= 4/3 - 7 - 3 - ye - «i + - -Ci - 3C2 , (42) 

and Eg is the gravitational binding energy of the body. This is known as the Nordtvedt effect 
(Nordtvedt 1968a), and can cause a difference in acceleration of the Earth and the Moon toward 
the Sun, and a resulting perturbation of the Earth-Moon orbit with a specific signature. Over 35 
years of lunar laser ranging have found no evidence for such an effect, and have placed the bound 
|??| < 9 X 10~^ (Williams et al. 2004). From the set of PPN parameter values for Whitehead's 
theory in Eq. (26), ?7whitehead = in strong disagreement with experiment. 



3.3. The binary pulsar 

Thirty years of timing of the binary pulsar 1913+16 have shown that its orbital period is 
decreasing at a rate Pb = —(2.4184 ± 0.0009) x 10~^^, in agreement with the GR prediction for 
gravitational radiation damping within a fraction of a percent (Weisberg &; Taylor 2005). Orbital 
damping has also been measured in the binary pulsars 1534+12 and the double pulsar 0737-3039AB, 
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again in agreement with GR. Unfortunately, Whitehead's theory has both the wrong sign - anti- 
damping instead of damping - and the wrong magnitude, ~ +4 x 10~^. The magnitude is so 
large because the reaction is a 1.5PN effect, rather that a f ^-timcs smaller 2.5PN effect, as in GR. 
One could change the sign of the effect, but not its magnitude, by assuming advanced, rather than 
retarded interactions. 

3.4. LAGEOS data 

Since 1992, precise laser tracking of two Earth-orbiting Laser Geodynamics Satellites (LAGEOS 
I and II) has made possible tests of general relativity in the vicinity of the Earth, in addition to 
its primary geophysical goals. Notably, the tracking data have been used to give a preliminary test 
of the "dragging of inertial frames", or Lense-Thirring effect, in which the rotating Earth causes a 
small precession of the planes of the orbits of the satellites. The NASA-Stanford Gravity Probe B 
experiment also aims to measure this effect with higher accuracy using orbiting superconducting 
gyroscopes. The effect depends on the PPN parameters 7 and ai, so both Whitehead's theory and 
GR agree on the prediction for this effect. However, the orbit of the LAGEOS II satellite has a 
small eccentricity, unlike LAGEOS I, and so its advance of perigee is also measured, along with the 
"nodal" precession of the orbit plane. 

Now, the multipole moments of the Earth's Newtonian gravity field also contribute to the 
nodal precessions and the perigee advance, indeed they overwhelm the relativistic effects. However, 
Ciufolini et al. (1997) found a particular linear combination of the three measurables, the two nodal 
precessions, $7/, and fi//, and the perigee precession of II, cj//, in which the effects of the leading 
I = 2 and I = 4 Newtonian multipoles would precisely cancel. The combinations depend on the 
known inclinations of the orbits relative to the equator. The uncertainties in the measured values 
of the remaining I > 6 multipoles then become part of the error budget of the experiment. 

The only difference in any of the relevant predictions between Whitehead and GR is the 
additional size-dependent term in the pericenter advance, Eq. (29). Because the LAGEOS II 

satellite is at two Earth radii, this can be a sizable effect (unlike the case with Mercury). Thus, 
the specific linear combination of predicted effects used by Ciufolini et al. gives the theoretical 
prediction (in milliarcseconds per year) 

^theory = + 0.29517/7 - 0.35^7/ 

= 60.2 - 109C + (errors) (43) 

where we have kept the PPN Whitehead parameter ^ but used the GR/ Whitehead values for 7, /? 
and ai, and where "error" denotes those due to the higher multipole moments. Using the actual 
tracking data, the measured value of this combination is A^xp = 66.6 milliarcseconds per year, plus 
measurement errors. The combination of all the errors leads to a total estimated error of about 25 
percent. Thus for the theory to match observation within 25 percent, the parameter ^ must lie in 
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the range 

-0.2 < ,e < 0.1 , (44) 

which thus excludes Whitehead's theory. It is hkely that this bound could be improved by maJiing 
use of dramatically improved Earth gravity models that have been derived from the GRACE and 
CHAMP geodesy space missions, which have reduced the errors in the Earth's multipole moments 
by significant amounts. 

3.5. Binary pulsars and momentum conservation 

The binary pulsar B 1913+ 16 provides an excellent system to test the momentum non-conserving 
effect described in Sec. 2.6, because it is highly relativistic, and because of the ability to do 
precise timing. For a moving system all measured periods will be offset via the Doppler effect 
(AP/P ~ v/c)\ accordingly, in an accelerating system periods will suffer a drift dP/dt ~ {a/c)P, 
and in a system with a changing acceleration, there will be a d?P/dt^ ~ {a/c)P. In the binary pul- 
sar, the center of mass acceleration predicted by Whitehead's theory changes because it is directed 
toward the periastron of the system, which rotates by 4 degrees per year. Indeed, in the 30 years 
since discovery, the center-of-mass motion (were it to exist) would have almost reversed itself. Yet 
precise timing of the pulsar 1913-1-16 has shown no evidence of any change in its spindown rate 
dP/dt, leading to an upper bound |d^P/dt^| < 8.5 x 10~^^ s~^ (Manchester et al. 2005). Using the 
neutron star masses and orbital elements inferred from the timing data, together with Eq. (4) of 
Will (1992), we find the predicted value d^P/df ~ 2.1 x 10^2^(2 cos CO s , where to is the periastron 
angle (we adopt the Whitehead value 03 = 0). With cos to varying between —1 and +0.59 over 
that period, we find the bound IC2I < 8 x 10~^. Notice that the mass values used were inferred 
using GR; in Whitehead's theory, it is conceivable that these values could be different from the GR 
values (as occurs in other theories that violate the strong equivalence principle). However, to evade 
this bound, either the inferred masses would have to be 10^ times smaller, or they would have to 
be the same to a part in 10®. This seems highly unlikely. 

4. Cosmological Considerations 

In addition to passing stringent tests at terrestial, solar system, and galactic scales, in order 
to be viable, a theory of spacetime and gravity must agree with the basic facts of cosmology: 
the expansion of the universe and the existence of the Cosmic Microwave Background. Of course 
neither was known when Whitehead formulated his theory. However at present, we are entering an 
era in which cosmological observations are becoming increasingly detailed and precise (Spergel et 
al. 2006). 

Already during the 1950's Synge (1954), using the spherically symmetric continuum version of 
Whitehead's theory developed by Rayner (1954), derived the form the Friedmann-Lemaitre metric 
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takes according to Whitehead. If r = ^—r}^i,x>^x^ ^ the curved metric is 

ds"^ = - - + + 7) ('^^^ + sinlix(d^^ + sin^ Odcf)'^)) , (45) 

with the density 

P=^, (46) 

and A = 8ttGK/9, with X a constant. Note that Synge's version of Whitehead's Universe, which 
has hyperbohc, fc = — 1 spatial cross sections, becomes empty and flat at late times, becoming more 
and more Milne-like. 

By contrast, current observations (Spergel et al. 2006) strongly indicate that our universe is 
currently of Priedmann-Lemaitre form with flat spatial sections and scale factor a(r) with jerk 
(Blandford 2004) 

d^a 



1, (47) 



dr^ 

and thus given by 

a(T) = sinh2/3[(3A/4)i/V] , (48) 



where A is the cosmological constant. As proper time r goes by, the universe is more and more 
accurately De-Sitter like, with 

a(r) = e(^/3)'/'- . (49) 



It seems that to be viable. Whitehead's theory requires, at the very least, a modification that 
incorporates the same effects as the cosmological term in Einstein's theory. The principal motivation 
behind Whitehead's alternative to Einstein's theory was the desire to retain fixed, non-dynamical, 
background-independent, causal relations between spacetime events which do not depend upon 
one's location in spacetime. Presumably, purely on the same aesthetic or philosophical grounds, 
one might argue that, as a fixed set of spatio-temporal relations, those of De-Sitter spacetime or of 
anti-De-Sitter spacetime are to be preferred to those of Minkowski spacetime since the underlying 
isometry groups in the former two cases are simple, rather than being a mere semi-direct product 
in the latter. Be that as it may, early on. Temple (1923) pointed out that this aim could just as 
readily be achieved by adopting the causal relations of a fixed maximally symmetric spacetime of 
constant curvature, e.g. a De-Sitter spacetime, as it could by insisting that they were the same as 
Minkowski spacetime. With this in mind. Temple sketched a generalisation of Whitehead's theory 
to incorporate a De-Sitter background which received enthusiastic support from Whitehead himself. 
An interesting Machian argument in its favour was made by Band (1929a), who pointed out that 
for positive cosmological constant it described a finite universe relative to which one could define an 
absolute acceleration. Actually Band claimed (Band 1929b) that Whitehead's theory was in gross 
violation of experiment. Later, Rayner (1955b) pointed out what he claimed were some errors in 
Temple's formulae. 
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Rather than recall the details of Temple's construction, which appears to have been almost 
completely forgotten, perhaps because the reference to it in Synge's influential reformulation of 
Whitehead's theory in modern notation (Synge 1952) is incorrect, we shall content ourselves with 
the remark that the obvious statement of the theory^ is that it amounts to linearising Einstein's 
theory with a cosmological constant around a De-Sitter background. This interpretation is consis- 
tent with Temple's finding that the perihelion advance agrees with that obtained by Eddington for 
the Schwarzscliild-Dc-Sittcr metric. If one accepts our interpretation, then the fact that the Kerr- 
De-Sitter solution is also of Kcrr-Schild form (Carter 1968) shows that Eddington and Rayner's 
observations (Eddington 1922; Rayner 1955a) may be extended to the full set of rotating solutions 
in a background De-Sitter spacetime. 

However, although incorporating a cosmological term may conceivably render Whitehead's 
theory in better accord with cosmological data, it will do nothing to alter the fact that it is in 
flagrant contradiction with observations at solar system and galactic scales, since the effects of any 
cosmological modification at these scales are negliglible. 
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